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ABSTRACT 

We determine the expected surface brightness and photometric signature of a white dwarf remnant 
population, issued from primordial low-mass stars formed at high redshifts, in today galactic halos. 
We examine the radial dependence of such a contribution as well as its redshift dependence. Such 
a halo diffuse radiation is below the detection limit of present large field ground-based surveys, but 
should be observable with the HST and with the future JWST project. Since the surface brightness 
does not depend on the distance, the integration of several galactic dark halos along the line of sight 
will raise appreciably the chances of detection. Both the detection or the non-detection of such a 
remnant diffuse radiation within relevant detection limits offer valuable information on the minimum 
mass for star formation in the early universe and on the evolution of the stellar initial mass function. 
Subject headings: Galaxy: halo - dark matter - galaxies: mass function - halo - evolution - stars: white 
dwarfs 



1. INTRODUCTION 

Modern theories of structure formation in the early 
universe suggest that galaxies have formed out of a pri- 
mordial generation of stars, with zero-like metal abun- 
dances, the so-called Population III stars. Our under- 
standing of the general properties of Pop III stars has 
improved significantly over the past recent years (see 
Barkana & Loeb 2001, Bromm & Larson 2004 for recent 
reviews). The mass distribution of this first generation(s) 
of stars, however, remains ill-determined. According to 
recent hydrodynamical models (Abel, Bryan & Norman 
2000, Bromm, Copi & Larson 1999, 2002) the first stars 
are very massive (>100M Q ), primarily because of the 
poor efficiency of the dominant cooling mechanism, due 
to molecular hydrogen (H 2 ) rotational transitions. Other 
calculations, however, find that the formation of Pop III 
stars following the very first-generation stars may extend 
to the ~ 1-10 M© range, depending on the protogalac- 
tic gas cloud conditions (Nakamura & Umemura 2001, 
2002), the effect of the intense radiation field (Omukai 
& Yoshii 2003) or on the efficiency of metal-enriched 
(Pop II. 5) star formation due to supernovae (Mackey, 
Bromm & Hernquist 2003, Salvaterra, Ferrara & Schnei- 
der 2003). In that case, primordial star formation may 
have led to an early generation of extremely or totally (in 
the absence of mixing of the ejected material in the SN- 
driven shell in the ISM) metal-depleted, long-live low- 
mass stars. The minimum mass for star formation in the 
early universe thus remains uncertain. This minimum 
mass, however, is likely to decrease with time, as (i) the 
virial temperature of massive collapsing halos will exceed 
T V i r > 10 4 , inducing radiative line-cooling from atomic 
hydrogen, (ii) the minimum ambient temperature, set up 
by the cosmic background radiation, will decrease and 
so will the thermal Jeans mass, and (iii) a small fraction 
of metals will be processed in the primordial short-lived 
stars, yielding a much more efficient cooling mechanism 
(Clarke & Bromm 2003, Bromm & Loeb 2003). When 
the transition between Pop 111+11.5, hereafter defined 
generically as our " primordial star generation" , to well 



identified Pop Il-like stars is occurring remains unknown, 
depending on ill-determined factors such as the effect of 
the ambient UV radiation on the main coolants or the 
efficiency of the medium heavy element enrichment due 
to the first SNII generations. 

The evolution of the initial mass function (IMF) from 
the (nearly-zero Z) Pop III to (low-Z) Pop II may in 
fact correspond to transition from a high-mass to a low- 
mass fragmentation mode of star formation, illustrat- 
ing the evolution of the IMF with structure formation 
in the universe. The recent discovery of a very metal- 
depleted ([Fe/H] = -5.3 ±0.2]) red giant (m w 0.8 M Q ) 
in our Galactic halo (Christlieb et al. 2002) indeed points 
to such a low-mass mode of primordial star generation. 
Moreover, the metal enrichment of clusters of galaxies 
suggests an IMF at high redshift with an average stellar 
mass larger than for present-day conditions, in the range 
~ 1-10 M-o (Portinari et al. 2004 and references therein). 
Such an evolution bears important consequences on the 
early history and the subsequent evolution of the uni- 
verse. While the UV radiation and a-element nucleosyn- 
thesis will be dominated by the more massive objects, 
a primordial, high-redshift stellar population extending 
to masses below about 10 M Q will (i) extend the nucle- 
osynthetic products to C, N and Fe enriched elements 
as well as s-processed elements and (ii) produce a pop- 
ulation of white dwarfs (WDs) sequestering part of the 
baryons. The quest for the direct detection of these rem- 
nants remains for now inconclusive, for present surveys 
do not reach deep enough magnitude or have a too small 
field of view to reach robust conclusions. An alternative 
observable signature of the presence of a primordial rem- 
nant population in today galactic halos would be their 
contribution to these halo surface brightness. N-body 
simulations (White & Springel 2000) suggest that the 
primordial star clusters might survive disruption during 
halo merging preceding galaxy formation and, through 
loss of kinetic energy by dynamical friction, end up in 
the central bulge, with low velocity. If, however, the 
low-mass stars are scattered by violent relaxation or tidal 
disruption of the cluster during halo merging, the rem- 
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nants of these early stars would be now incorporated into 
galactic external halos. Which of these two processes 
dominates remains to be determined unambiguously and 
depends on various factors such as the mass ratio and 
the compactness of the mergers or the efficiency of the 
SN-driven winds in the clusters. In the present paper 
we will assume that the primordial remnant population 
follows the isothermal density distribution of the dark 
halo. If these remnants are located in the central bulge 
or spheroid, their surface brightness will be largely dom- 
inated by the main sequence (M-dwarfs) contribution. 
In that case only spectroscopic and high-proper motion 
detections of individual objects will allow their identi- 
fication as PopIII remnants. The nature of these pri- 
mordial remnants (black holes, neutron stars or white 
dwarfs) carries the signature of the minimum mass of 
the IMF at the early epoch of star formation in the uni- 
verse. Since the surface brightness does not depend on 
the distance (see below), the integration along the line of 
sight of several galactic dark halos adds up each contri- 
bution, yielding eventually a detectable signature. Both 
the observation or the non observation of such a stellar 
relic radiation would bring important information on the 
low-mass end of the stellar mass spectrum and thus on 
the dominant process of star formation in the early uni- 
verse. This is the aim of the present paper to examine 
this issue. 

2. VARIATION OF THE INITIAL MASS FUNCTION 

As mentioned above, several arguments suggest that 
the IMF was more top-heavy at high redshift than in 
present day environment. Figure ^ illustrates the pos- 
sible variation of the IMF from PopIII+II.5, formed at 
high-redshift (z>10 or so) with very low metal enrich- 
ment ( [M/H] < -2), to PopII stars, with [M/H] 2.0 

to -1.0, characteristic of the spheroid and globular cluster 
population, to solar metallicity PopI stars representative 
of the disk population. These IMFs correspond to eqns. 
(17), (20) and (24) (see also Tables 1 and 2) of Chabrier 
(2003), respectively. These IMFs are best described by 
a power-law at large masses rolling down to a lognor- 
mal form below some characteristic mass, a behaviour 
expected from turbulence-driven fragmentation (Padoan 
&: Nordlund 2002), but can be approximated by the fol- 
lowing form, as initially suggested by Larson (1986) : 



(m) am" exp 



,moyj 
m 



(1) 



As mentioned above, it is unclear whether a WD- 
dominated primordial IMF such as the one displayed in 
Figure H is representative of the primordial PopIII+II.5 
stellar generation. However, as argued above, different 
calculations suggest that the next to very first stellar 
generation may extend into this domain. Moreover, the 
aforementioned recent observation of a 0.8 M Q red gi- 
ant with a strong iron depletion ([Fe/H] — —5.3) sug- 
gests indeed that the primordial IMF quickly evolved into 
the WD progenitor mass range. In order (i) to maxi- 
mize the contribution from the putative WD-progenitor 
population and thus estimate the maximum constraint 
on this population, and (ii) to examine the dependence 
of the diffuse radiation upon the IMF and thus WD 
mean mass (jtiwd), we chose for illustration the IMF 




Fig. 1. — Illustration of the possible variation of the IMF with 
redshift and characteristic stellar populations. The PopI, PopII 
and PoPIII+II.5 IMFs correspond to Table 1, Table 2 and eqn.(24) 
of Chabrier (2003), respectively. The dotted line corresponds to the 
Salpeter IMF. All IMFs normalized to °° m (dN/dm) dm = l. 



suggested by Chabrier, Segretain & Mera (1996), with 
mo = 2.0 M , a = 5.1, (3 = 2.2 (their IMF1), yielding 
( TOwD ) = 0.62 M and m = 2.7M , a = 5.7, = 2.2 
(their IMF2), which yields (m WD ) = 0.68 M Q . 

3. EXTRAGALACTIC BACKGROUND LIGHT 

An important constraint on the presence of such a pri- 
mordial remnant population is the contribution of its 
progenitors to the energy content of the universe, i.e. 
to the cosmic background light. This contribution was 
first estimated by Madau and Pozzetti (2000). These 
authors, however, used solar-metallicity or only slightly 
metal-depleted (1% Z Q ) stellar evolution models to eval- 
uate this contribution. For a given mass, zero-metallicity 
stars are brighter and thus put more severe constraints 
on the related energy background. More recently San- 
tos, Bromm & Kamionkowski (2002) carried out more 
detailed calculations. They considered, however, only 
very massive stars (m>300M Q ) and found out that the 
spectral properties of these stars are essentially indepen- 
dent of mass, when normalized to unit stellar mass, so 
that the integrated spectrum of such a population de- 
pends only on its total mass, not on the IMF. This re- 
sult does not hold anymore for lower mass stars, and has 
been shown recently to be partially incorrect even for the 
highest mass stars (Schaerrer 2002). More recently, Sal- 
vaterra & Ferrara (2003) and Magliocchetti, Salvaterra 
& Ferrara (2003) have shown that Pop III stars formed at 
z > 9 can contribute appreciably, possibly entirely to the 
detected residual (i.e. after galaxy contribution subtrac- 
tion) NIRB radiation, without major conflict with the 
IGM observed metallicity limits. All these calculations 
suggest that a substantial fraction of the collapsed (viri- 
alized) baryons, ~ 10 to 60% depending on the primor- 
dial IMF and the redshift range of formation, may have 
been turned into this primordial generation of stars. 

This point is illustrated below. Following Madau & 
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Pozzetti (2000), we calculate the bolometric extragalac- 
tic background light per steradian Iebl observed at 
Earth today (z = 0), due to a burst of formation of pri- 
mordial stars of zero-metallicity at time tp, i.e. redshift 
zp. Since zero-metallicity stars are brighter and hot- 
ter than their metal-enriched counterparts, they evolve 
more quickly and the redshift dependence of the pro- 
duced light is likely to differ from the estimate of Madau 
& Pozzetti (2000). A star-burst approximation is rea- 
sonable since primordial star formation is believed to 
have occurred somewhere between zp ~ 20-30, when 
small scale halos were formed from primordial fluctua- 
tions, and Zp ~ 6, about the end of the reionization. In 
the standard cosmological model adopted in the present 
calculations, (/i, f2 m , Qa) = (0.70,0.30,0.70), this spans 
~ 10 8 yrs, a fairly short time compared with the age of 
the universe today, th = 13.7 Gyr. The Iebl is given by: 

_ c /*» Pbol (t) 
i-EBL - -j— / — — j-rdt (2) 
4?r J tF l + z{t) 

where 

Pboi{t)= f L(r)p(t-T)dT (3) 
Jo 

is the total bolometric luminosity density for a stellar 
formation rate per comoving volume p and 

L(t) = / l(m, T)<p(m)dm (4) 

J mint 

is the luminosity produced by a stellar population of 
same age r, from a minimum mass m,i n f to a maximum 
mass m sup . We use the zero-metallicity evolutionary se- 
quences of Marigo et al. (2001), with m, n / = 0.1 M Q and 
nisup — 1000 Mq. The term l(m,t) denotes the luminos- 
ity of a star with mass m at age r, and <j)(m) — dN/dm 
is the IMF normalized such that the total stellar mass 
M = f™ sup m(f>(m)dm = 1. For a burst SFR, eqn.© re- 
duces to pboi(t) — p mt L(t — tp), where p mt is the mass 
density under the form of stars. 

Figure displays Iebl{zf) expected for the afore- 
described population of primordial stars as a function of 
redshift of formation zp for various mass fractions X+(z) 
of baryons under the form of stars at redshift z : 

n*(z) = = X*(z) x fi B , (5) 

Pc 

where pc = 3Hq/8ttG is the universe critical density 
and Qb is the baryon density in units of pc- The ob- 
served (frequency-integrated) brightness of the cosmic 
background radiation today is Iebl — 100 nW m~ 2 sr _1 , 
with a lower limit Iebl — 45 nW m~ 2 sr _1 (Hauser & 
Dwek 2001), indicated by the dotted line in the figure. 

These calculations are by no means as accurate as 
the detailed calculations mentioned earlier. In partic- 
ular they do not incorporate other contributions like e.g. 
normal galaxies. However, they illustrate the fact that a 
significant mass fraction of primordial stars, possibly as 
much as fl+ = 0.5 x fi^ can be formed at redshift z > 5 
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Fig. 2. — Extragalactic background light received at Earth today 
due to a burst of zero-metallicity stars formed at redshift zp with 
IMF1 and containing a mass fraction X+ of the nucleosynthctic 
baryon density. The dotted line indicates the lower limit of the 
observed light. 



and contribute a diffuse background light today well be- 
low the observed lower limit. Since the IMF extends in 
the WD-progenitor range and since WD progenitors re- 
turn about 70% of their mass to the ISM with the IMF 
Q (assuming a similar WD-mass vs progenitor-mass re- 
lationship as for the disk WDs), this implies that, in prin- 
ciple, as much as ~ 15% of the nucleosynthetic baryon 
density can be trapped in a population of remnants from 
a primordial generation of low-mass stars. 

This is of course a large overestimate of the true frac- 
tion of baryons possibly processed into PopIII+II.5 rem- 
nants. This fraction can be estimated as follows. The 
fraction of collapsed baryons at z ~ 10 is Fb ~ 0.1 
(Haiman & Loeb 1997 1 , Abel et al. 1998, Bromm ct 
al. 1999). The star formation efficiency in today star- 
forming regions range from /* ~ 0.1 to 0.3 and, as 
mentioned above, may reach ~ 0.6 at high red- 
shift (Salvaterra & Ferrara (2003), Magliocchetti et al. 
2003). Given an estimate of f esc ~ 0.1 from the obser- 
vation of local galaxies, where f esc is the escape frac- 
tion of ionizing photons within these protogalaxies, this 
is consistent with the constraint /* x / esc <0.1 from 
the WMAP Thomson optical depth and the lack of H 
I Gunn-Peterson through at z <; 6 (Fig. 5 of Wyithe 
& Loeb 2003). This yields a fraction of BBN bayons 
X* = /* X Fb(z) w 1-10% possibly processed into 
PopIII+II.5 stars at large redshift. Recent studies sug- 
gest that the direct detection of the radiation of these pri- 
mordial star clusters, responsible for the universe reion- 
ization, might be in reach of the JWST or even the HST 
missions (Venkatesan & Truran 2003, Stiavelli, Fall & 
Panagia 2004). On the other hand, with the type of IMF 
mentioned earlier, about 30% of these baryons might be 
under the form of WD remnants in galactic halos today. 

1 Haiman & Loeb (1997) simulations are based on Press- 
Schechter SCDM and different cosmological parameters but mod- 
ern simulations do not seem to change drastically these results. 



4 



The glow of primordial remnants 



The photometric signature of such a remnant population 
in galactic halos is examined below. 

4. DIFFUSE EMISSION OF POPULATION III WHITE 
DWARFS 

The mass fraction Xwd = Pwd/pHqi where ph e — 
9 x 10~ 3 M© pc -3 is the local dark matter density, of WD 
remnants in the Milky Way halo, issued from the primor- 
dial low-mass star generation considered in the previous 
section, remains presently undetermined. The value in- 
ferred from microlensing experiments has been steadily 
decreasing over the years, and is now estimated to be be- 
tween X WD = and —20% (Alcock et al. 2000, Afonso 
et al. 2003), with X WD < 5% at the 95% C.L. from 
the EROS team. A Xwd ~ 2% population of these halo 
WDs has been claimed to have been identified recently 
from kinematic and spectroscopic observations (Oppen- 
heimer et al. 2001). Although this result remains highly 
controversial (see e.g. Reid, Sahu & Hawley 2001, Reyle, 
Robin & Creze 2001, Flynn, Holopainen & Holmberg 
2003), a detailed analysis of the age and kinematic prop- 
erties of these objects seems to confirm that part of the 
sample does belong to the Galactic halo (Hansen 2001). 
Present surveys, however, either have too small field of 
views or do not reach faint enough magnitudes to detect 
the peak of the halo WD luminosity function and thus the 
majority of such a halo WD population remains presently 
out of reach (Chabrier 1999a Figure 1). As pointed out 
by Chabrier (1999a), the only robust constraint yielding 
an upper limit for Xw jy (besides the maximum baryonic 
fraction estimated in §3), arises from stellar nucleosyn- 
thesis of helium and heavy element production. Recent 
detailed calculations seem to exclude a value Xwd > 1% 
to 10% (Fields, Freese, & Graff 2000, Brook, Kawata & 
Gibson 2001). In any events, a primordial, high- velocity 
(u_L>100 km s _1 ) WD population with Xwd > 1%> as 
suggested by Oppenheimer et al. (2001), would exceed 
significantly the WD population expected in the Galac- 
tic spheroid, pwd — (1.8 ± 0.8) x 10~ 5 M Q pc~ 3 , i.e. 
X W D aph ~ 0.2 ± 0.1% (Gould, Flynn & Bahcall 1998, 
Chabrier 2003 Table 3), and would thus imply a primor- 
dial IMF different from the one characteristic of today 
environments. At any rate, the identification of any halo 
WD population would bring precious information about 
the minimum mass for star formation at early stages of 
galactic formation. It is thus interesting to explore the 
expected diffuse emission of such a population. 

Following the formulation of Adams & Walker (1990), 
we calculate the radiation signature of a WD population 
with mass fraction Xwd in an isothermal halo. The 
specific luminosity per unit mass (r„) A in a broadband 
filter characterized by central wavelength A and width 
AA of a WD population obeying the distribution given 
by an IMF dn/dmwD, reads: 



and 
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Here, (mi n f,m sup ) = (0.5 M , 1.2 M Q ) denote respec- 
tively the minimum and maximum masses for a stan- 
dard WD mass spectrum, whereas L Vx {mwD,T~) — 
4:TrR WD (mwD,'r)(J r v)\('mwD,'r) is the luminosity per 
unit frequency averaged over the filter A for a WD of 
mass mwD and age r, emitting at its surface a flux T v . 

The surface brightness l v and the flux density at fre- 
quency v per unit mass f u — T v /A-ks 2 observed at Earth 
from such a population emitting at distance s in a galac- 
tic halo of density profile p[r(s)], where r and s denote 
distances from the Galactic centre and from the Sun, re- 
spectively, reads: 



h = 



p[s(b,l)]dV 



f u _ Xwd f r„ 



= x wd I p[s(b,l)]ds 



(9) 



The apparent and absolute magnitudes in the filter A 
read, respectively : 

m A = -2.5 1og(/ !/ ) A + C A (10) 
A/ A = -2.5 1ogZ„ A + 2.5 log(47r[10pc] 2 ) + C A (11) 

where (f v ) x = J °° fvS\(i/) di>/ J °° S\(v) dv is the flux 
distribution in the broadband filter characterized by the 
transmission function Sx^), and C\ denotes the nor- 
malization constant plus zero-point correction in the fil- 
ter. The present calculations are done in the stan- 
dard filters, with magnitudes normalized to Vega (Bessell 
& Brett 1988, Allen 1991). The absolute magnitudes 
M\(mwD,T~) and radii RwD(mwD, t) of the WDs are 
taken from the cooling sequences of Chabrier et al. 
(2000) 2 , based on the spectral energy distributions of 
Saumon & Jacobson (1999). This describes the evolution 
of so-called DA WDs, with hydrogen rich atmospheres, 
the most likely population, given WD accretion rates 
(Chabrier 1999b). At any rate, DA WDs, for a given 
age, will be order of magnitude brighter than helium-rich 
atmosphere so-called DB WDs, and will thus dominate 
the diffuse radiation flux. Note that the age r of the WD 
includes both the WD cooling time t coo i and the main se- 
quence lifetime tMS of the progenitor of mass m, with the 
condition that the WD population today, at halo age th , 
must obey the condition (Chabrier et al. 1996, Chabrier 
1999a): 



th = tMs(m) + t coo i(L, vtlwd) 

4.1. Milky Way halo 
The surface brightness 10 is rewritten: 



r 

l v = Xwd -r-Ph @ R®Z\pi 



(12) 



(13) 



WD 



Available at http:/ /perso. ens-lyon.fr/gilles.chabrier/WD2000. tar 
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where Rq = 8.5 kpc is the Galactocentric position of the 
Sun, and l(b, I) = (p/jq-Rq)" 1 ^ mal p[s(b,l)]ds denotes 
the angular integration along the line of sight. 

The density distribution of a flattened isothermal halo 
with galactocentric cylindrical coordinates (R, z) reads : 



R% 



P{R,Z) Ph @ (- qAlccQsq ; R 2 + a 2 + z 2/ q 2 



4 7r G q Arccos q R 2 



z 2 / q 2 



(14) 



where a ~ 5 kpc denotes the halo core radius and q ~ 0.7 
the oblateness. 

For such a density distribution, the integral 1(13(1 can 
be evaluated analytically: 



Table 1 

Magnitude (in mag/arcmin 2 ), top row, and surface 

BRIGHTNESS I v (iN mJy/STR), BOTTOM ROW, FOR A 
REMNANT WD POPULATION WITH X wd = 1% IN THE MlLKY 

Way halo at (6,11) = (45°, 0°), for a Dirac IMF. The 

FLUX AND MAGNITUDES ARE CENTERED AT Aq AND 
AVERAGED OVER AA (iN fjM) . 



Filter 


B 


V 


R 


I 


J 


H 


K 


M 


Ao H 


0.44 


0.55 


0.70 


0.90 


1.25 


1.65 


2.20 


5.0 


AA [fjm] 


0.10 


0.09 


0.22 


0.24 


0.30 


0.35 


0.40 


0.3 


m wd = 0.6 M Q 


28.7 


27.6 


27 


26.4 


25.9 


26 


26 


25.6 




180 


374 


531 


742 


890 


307 


280 


121 


m wd = 0.8 M Q 


29.6 


28.5 


27.8 


27.2 


27 


27.1 


27.2 


26.7 




81 


171 


264 


351 


343 


107 


93 


44 



1(6,0 = 



-{tan-(f) 



- U p-Y max (l + Q 2 sin 2 b) 
— tan ( 



)} (15) 



with fj, — cos b cos I 
VW - 1, and £ - [(a 



a = a/i?Q, 
2 < 1)(1 ' 



Y = s/R G 
Q 2 sin 2 b) - /i 2 ] 



Q = 

r- 2 ] 1/2 - 

Equation 1)15(1 recovers equation (9b) of Adams & Walker 
(1990) for a spherical halo (q = 1). 

The maximum extension R max of the halo is deter- 
mined by the total mass Mh and the asymptotic velocity 
v rot of the considered halo: 



Rr, 



G^(^ 



q Arccos q 



(16) 



where v rot « 220 km s^ 1 and M^2x 10 12 M Q for the 
Milky Way. 

A first estimate of the diffuse emission can be ob- 
tained by assuming a Dirac function for the IMF: 
dN/d(m,TNT>) — N Q x 5(m — m Q ), where N Q is fixed by 
the normalization condition N = Xwd x (-M/,/mo). In 
that case the specific luminosity simply reads : 



T v = 



L v {m ,' 
m 



(17) 



Table I displays the expected surface brightness and 
apparent magnitudes of the WD diffuse emission in the 
Milky Way dark halo, at (b,l) = (45°, 0°), for such a 
Dirac IMF with m = 0.6 M Q and m = 0.8 M Q , re- 
spectively, for Xwd — 1% and t/, = 13 Gyr. Figure 
illustrates the apparent magnitudes of the WD dif- 
fuse emission for the same conditions for different IMFs, 
in various optical and near-infrared bands. The mass 
functions IMF1 and IMF2 refer to the WD-dominated 
IMFs of Chabrier et al. (1996) (eqn.Q}) mentioned in 
§2. The IMF PopIII+II.5 displayed in Figure 1 gives re- 
sults similar to IMF2. The IMF labeled "disk" refers to 
the IMF of the present Galactic disk (Chabrier 2003, Ta- 
ble I). The dependence of the flux upon the IMF between 
IMF1, which yields (m WD ) ~ 0.6 M , and IMF2, which 
yields (m W D) — 0.7 M Q , illustrates the strong depen- 
dence of WD cooling upon the WD mass (see Chabrier et 
al. 2000), whereas the difference between IMF1 and IMF 




Fig. 3. — Surface brightness of the Milky Way halo containing 
a WD fraction X\yd = 1% of the dark matter density in different 
passband filters for three different IMFs, as described in the text. 



"disk" , which both have similar WD average mass, stems 
from the larger average WD luminosity (L v ) in the case 
of IMF 1. Note from equation ((15(1 that the received radi- 
ation has an angular dependence: magnitudes are found 
to increase by about 1.5 mag from b = 0° to b > 90° and 
by about 1 mag from I = 0° to I = 180°. As mentioned 
above, the magnitudes have been calculated in standard 
filters, normalized to Vega. The flux of Vega decreases 
monotonically from ~ 0.5 /im (V-band) redward, yielding 
as a spurious effect decreasing WD magnitudes redward 
of the J band. The flux, however, peaks in the I and J 
bands, as seen in Table I. These bands thus seem to be 
favored for the observation of the remnant WD surface 
brightness. Since I v is proportional to Xwd (eqn.JHJl), 
these curves can be easily scaled to different values of the 
expected halo WD mass fraction. Although difficult to 
detect with present instruments, the detection of such a 
halo WD emission might be observable with the HST or 
with the JWST future satellite mission. We will come 
back to this point later. An other possibility, as explored 
in the next section, is the observation of halos of external 
galaxies. 
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4.2. External galaxy halo 

As seen from eqn.@, the surface brightness does not 
depend on the distance, because of the canceling r 2 de- 
pendence of the flux density and volume integration. We 
thus calculate now the expected surface brightness of our 
primordial remnant population for the halo of an external 
galaxy, located at distance d from Earth. The galactic 
density distribution is supposed to obey the same rela- 
tion as given by ean. (|14|l : 

P^ z ^P'^ R , + J +z2/q , (18) 

We assume for this galaxy a rotation curve, and thus 
a dynamical halo identical to the Milky Way conditions, 
with a similar normalization p' h — ph @ x (Rq + a' ) , 
a' = a, q' = q. Such conditions apply for example to the 
spiral galaxy NGC5907 (Lequeux et al. 1998, Figure 2). 

We note x, y, z the cartesien galactocentric coordinates 
of the external galaxy, with z the direction perpendicular 
to the plane of the galactic disk. Galaxies seen edge-on 
are more favorable targets since in that case observations 
far away from the galactic center are not polluted by the 
light from the disk. In that case, the surface brightness 
I' v observed at Earth now reads, with x denoting the 
distance along the Earth-galaxy line of sight and V Q b s 
the observed galactic volume: 

I » = n=Uv obs dVp{x,y,z,) ^ 

c +x ' r 
= L dx{d + xfp[x -' mz ' ) MdfxT 

= ^P' h ^Uz), (19) 

with 



Table 2 

Same as Table 1 for the cumulative contributions 
of the Milky Way halo and a similar external halo 
(namely NGC5907) at 40 kpc from the galactic 
center, for the IMF1 and Xwd = 1%. At this 

DISTANCE, THE MlLKY WAY HALO CONTRIBUTES ~ 70% TO 
THE TOTAL FLUX. 



B 


V 


R 


I 


J 


H 


K 


M 


mag 28.2 


27.4 


26.9 


26.4 


26.0 


25.9 


26.0 


26.2 


I v (mjy/str) 284 


454 


570 


732 


831 


314 


302 


71 



24 



80 
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Fig. 4. — Surface brightness radial profile from galactic center 
in different bands for an external galactic halo composed of a mass 
fraction Xwd = 1% OI primordial WDs, including (solid line) and 
without (dotted line) the Milky Way halo contribution. 



dx 



2 



yV + z 2 /q' 2 + a' 2 



x 2 + y 2 + z 2 /q' 2 + a' 2 

tg- 1 



\/y 2 + z 2 /q' 2 + a' 2 



(20) 

and the maximum halo extension R' max is fixed by the 
halo mass, as in eqn.(Hni). For M h — 2 x 10 12 M Q , 
R' max ~ 150 kpc, so that for galaxies situated at dis- 
tances d ~10-100 Mpc, the halos extend to angular sizes 
~ 0.1°-1°. Figure 21 portrays the expected radial de- 
pendence in the galactic plane of the halo of the galaxy 
NGC5907 alone ((b, I) = (51.1°, 91.6°)) and of the cumu- 
lative contributions of this halo plus our own halo, for 
X w d = 1%, in various passbands. The total flux values 
are given in Table II. Depending on the radial distance y, 
our own halo contributes from ~ 20% near the galactic 
center to ~ 80% at 80 kpc to the total flux. The diffuse 
emission is below the detection limit of present large, 
deep-field cameras like MEGACAM or WIRCAM, but 
should be in reach of space-based HST or JWST deep 
field surveys. Note that both the angular modulation 
and radial dependence of the signal will help separate 



the halo radiation from other possible sources, helping 
correcting for the sky level contribution. Present obser- 
vations of NGC5907 (Lequeux et al. 1996, Rudy et al. 
1997, Zheng et al. 1999) extend only to 100" from the 
galactic center, i.e. ~ 7 kpc for a distance d = 14 Mpc. 
At this distance, contamination from the disk and the 
bulge is not negligible. 

Alternatively, an indication on the primordial WD 
population of dark halos would be their contribution at 
high redshift. High redshift means younger ages and thus 
brighter WD luminosities. Figure [3] displays the redshift 
dependence of the halo WD I-band surface brightness, for 
various redshift values up to z = 5, i.e. th ~ 1 Gyr, for 
Xwd — 1%- Already at z = 1, th ~ 6 Gyr, the remnant 
emission is about 1 magnitude brighter than at z = 0. 
At high redshift, however, some of the WD progenitors 
are still on the main sequence or the AGB phase and 
contribute significantly to the background light. A fully 
consistent calculation requires chemico-evolution codes 
(see e.g. Lee et al. 2004). Appropriate color-cuts, how- 
ever, might be able to observe this contribution from 
halo WDs at larger redshifts. Figure illustrates the 
color dependence of the WD halo signature under the 
same conditions. 
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Fig. 5. — I-band surface brightness radial profile for an external 
galactic halo with a fraction Xwd = 1% of primordial WDs for 
different rcdshift values. 




2 3 
z [redshift] 

Fig. 6. — Color dependence as a function of rcdshift for a galactic 
halo composed of a mass fraction Xwd = 1% of primordial WDs. 



5. CONCLUSION 

In this paper, we have examined the possibility for a 
fraction of the baryons in the universe to be trapped 
in the WD remnant population of a primordial genera- 
tion of intermediate or low-mass stars, formed before the 
end of the reionization period, i.e. in the redshift range 
5 < z < 20. For such redshift values, a non-negligible frac- 
tion of the nucleosynthetic baryons can be in the form 
of such WD progenitors without producing unacceptable 
light today. As a possible signature of such relics, we 
calculate their expected surface brightness both in our 
own halo and in external galactic halos. The expected 
radiation lies below the detection limits of ground-based 
present surveys, but the detection of this glow should 
be in reach of the HST or the JWST, in particular if 
several external galactic halos can be observed along a 
given line-of-sight. Both a successful or a null detection 
of such a diffuse light, within relevant detection limits, 
will provide precious information about the mode of star 
formation and the evolution of the IMF at early stages 
of galaxy formation. The background light contribution 
of these WDs is found to increase appreciably with red- 
shift, suggesting possible detection in dedicated surveys, 
with appropriate color cuts. Contributions from red gi- 
ants and main sequence stars at these redshifts, however, 
must be calculated consistently with chemico-evolution 
models, in order to examine the favored detection bands. 
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